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We investigate the deconfinement transition in SU(iV c ) gauge theories, and properties of the de- 
confined phase. A detailed lattice study of SU(4) and SU(6) gauge theories are conducted, and 
finite volume and cutoff effects on thermodynamic observables are studied. The scaling of the 
deconfinement transition point with lattice spacing is used to calculate A^j. The continuum es- 
timates of the thermodynamic quantities are used to study properties of the gluon plasma. In 
particular, the approach to conformal limit is studied. We do not find any evidence of a strongly 
coupled, conformal phase in these theories. 
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1. Introduction 

At very high temperatures, strongly interacting matter is known to exist in a deconfined, chi- 
rally symmetric state. The nature of the transition to this state is quite sensitive to the quark sector. 
For physical quarks, one has no phase transition, only a crossover. For very light quarks one ex- 
pects a chiral symmetry restoring transition, whose order will depend on the mass and the number 
of quark flavors. On the other hand, for pure SU(3) gauge theory, which is the limit of QCD 
with infinitely heavy quarks, one has a first order, deconfming transition. See Ref. [jl]| for further 
discussions of the phase diagram. 

While the nature of the transition is sensitive to the quark content, one can expect similarities 
in the properties of the high temperature phase of the pure gauge theory and of QCD. In particular, 
while at very high temperatures, the theories are expected to be weakly coupled due to asymptotic 
freedom, at moderately high temperatures ~ 1 — 3T C , large deviations from the weak coupling 
theory are seen in both theories. An understanding of the interactions in the pure gauge theory 
plasma will surely help in our understanding of the QCD plasma. 

A study of the SU(iV c ) gauge theory in the limit of large N c may be of help in this respect. Large 
N c arguments have been used to understand various features of the SU(3) deconfinement transition 
[§]. Also, interesting phase structures have been suggested for SU(iV c ) theory with quarks, when 
the baryon number chemical potential }Xb ^ N c Analogies have also been drawn between 
the high temperature phase of QCD and the analytically tractable, JV =4 supersymmetric SU(°°) 
theory Of course, connecting the S\J(N C — > °°) theory to SU(3) will require that the corrections 
in 1 /N c are small. An estimate of the size of the corrections can be obtained by studying SU(iV c ) 
theories with different N c on lattice. 

In this paper, we report on a study of the S\J(N C ) gauge theories at finite temperatures for 
N c = 4,6. In Sec. ^ we study the deconfinement transition. Using a finite volume analysis at 
different lattice spacings, we establish the first order nature of the transition. We also study whether 
asymptotic scaling holds in the transition regime, and use it for setting the scale through T c /A^. 
More details about the material in this section can be found in Ref. In Sec. [3] we study the 
equation of state of the SU(jV c ) gluon plasma. A comparison of the energy density and pressure for 
SU(3,4,6) gives us an idea of the size of the 1 /N c corrections. In the final section we discuss the 
physical implications of our results. Some of these results were earlier reported in ref. |jlj. 

SU(Af c ) gauge theories have been investigated previously on the lattice. In particular, the de- 
confinement transition has been studied in Ref. ^ [7j], while the equation of state close to T c has 
been studied in Ref. Our work focuses on the continuum limit, and also explores higher 

temperatures in the equation of state, which allows us to investigate the issue of approach to con- 
formality. Equation of state at high temperatures have also been reported in this conference by 
Panero [0], who used a coarser lattice, but reported also results for SU(8). 

2. T c and scale setting 

SU(3) gauge theory is known to have a weak first order transition, with the Polyakov loop, 

1 1 N > 
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acting as the order parameter. Here V = is the spatial volume and is the temporal extent of 
the finite temperature lattice, and U^(x,xo) is the gauge link matrix connecting the sites x = (x,xo) 
and x + afX, where a is the lattice spacing. 

We have carried out a detailed study of the deconfinement transition in SU(4) and SU(6) gauge 
theories, on lattices with spacings of a -1 < 6T C with the Wilson action. Finite volume analysis was 
carried out on lattices with N t = 6 and 8 by using various aspect ratios, £ = N s /N t , in the range ~ 
2.5 - 4. In the transition regime, L shows a clear signal of a discontinuous transition, tunnelling 
between the low temperature phase around L = and N c non-zero values, related by rotations of 
2n/N c , for the high temperature phase. Figure |] shows the histogram for \L\ for SU(4) gauge theory 
on N t = 8 lattices, showing clearly the tunnelling between the confined and deconfined phases, with 
the two-peak structure getting sharper as one goes to larger lattices. This is the expected behavior 
for a first order transition. For a more quantitative analysis, we look at the susceptibility of the 
order parameter, 



%L = Nl{(\L\ 2 )-(\L\) 2 } 



(2.2) 



For a first order transition, %l ~ V. In Fig. [I] we also show the result of a multihistogramming 
analysis of Xl/V for the same set, confirming the first order nature of the transition. More details, 
and figures for SU(6), can be obtained in Ref. [B]. 
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Figure 1: Finite size analysis of the deconfinement transition for SU(4) gauge theories, on Nt = 8 lattices. 
(Left) Histogram of \L\, for lattices with N s = 22- 30. (Right) Volume dependence of the |L| susceptibility. 



To confirm that one is looking at a physical transition, one needs to confirm that /3 e (A^) has 
the proper scaling behavior, so that the transition temperature, T c = l/N t a(fi c ), is independent of 
N t . To study this scaling, we also add runs at N t = 10 for both SU(4) and SU(6), and N t = 12 for 
SU(4), at a single volume (N s = 24). We use the two-loop renormalization group evolution (RGE), 
and use j3 c (A^) for the different N t to set a temperature scale. In Fig. || we show the temperature 
scale for N t = 6, set by using the p c (N,) for different N t , for both SU(4) and SU(6), Here we have 
also included j8 c (A^ = 5), taken from Ref. As the figure reveals, while some deviations from 
scaling are seen at smaller N t , in particular for ,/V, = 5, very good scaling is observed for N t > = &. 

This allows us to use the data for N t > 8 to get an estimate for the scale, Ayg, of the theory [ffO]. 
We use the known expansions of plaquette to get the renormalized coupling, and use the two-loop 
running to convert it to A^. Since we are performing the calculation at a finite order of the weak 
coupling expansion, we also investigate the scheme dependence, by calculating the renormalized 
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Figure 2: Temperature scale, T/T c (fi), for N, = 6, using p c (N t ) for different JV,, for (a) SU(4) and (b) SU(6). 



N c 


E-scheme 


V-scheme 


MS-scheme 


3 


1.19(3) 


1.12(3) 


1.20(2) 


4 


1.235(1) 


1.153(1) 


1.236(1) 


6 


1.222(1) 


1.135(1) 


1.217(1) 



Table 1: Continuum estimate of T c /A^ for SU(N C ) gauge theories with N c = 3,4 and 6. The renormalized 
couplings in different schemes are obtained from plaquette data, and converted to aA using the two-loop 
RGE. These are converted to T c /A^ using the known relations between schemes. A constant fit to data for 
N t > 8 is used for the continuum values. 



coupling in some other scheme, getting the A parameter in this scheme, and then converting it to 
Ajjj. In table [T] we show the results of such estimates, where we have also used the E-scheme 
[11] and the V-scheme [p^J. The data for N c = 3 are taken from Ref. [|l3|]. We find that there is 
considerable scheme dependence in the estimation of T c /A^. We also see that within a scheme, 
the dependence on the number of colors is quite weak. 



3. Equation of State 



For estimating the thermodynamic quantities, we use the integral method P13fj. The pressure, 
p, and £ — 3/7, where e is the energy density, are calculated from the plaquette data, [13] 

'ft 



3p 



f4 



6N?^AP(P,T) 
da 



(3.1) 



where AP(fi,T) is the difference in the plaquette observables between the finte temperature lat- 
tice and the corresponding zero temperature (symmetric) lattice, calculated at the coupling /3, 
AP(P,T) = P(f},T)-P(J},T = 0). T Q is some reference temperature. We find that p(T) ~ 
within our errors till temperatures ~ 0.9T C , and evaluate p/T by taking fio(T < O.ST c ) as the 
lower limit of the integral. For we use the two-loop result, with the coupling defined through 
the V-scheme. 

The results of Sec. || indicate that two-loop beta function should work fine when a < l/ST c , 
while some cutoff effect may be seen in coarser lattices. Since we are interested in the thermody- 
namic and continuum limit results, we investigate the cutoff and finite volume effects by using two 
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Figure 3: Volume and cutoff dependence on (e - 3p) /T 4 (left) and p /T 4 (right) for SU(4) (top) and SU(6) 
(bottom) gauge theories. The results shown have been normalized by the known discretization error in free 
theory, so that the Stefan-Boltzmann limit for p/T 4 for both the N t is — l)7T 2 /45 (see text). 



aspect ratios and two cutoffs. Note that the Stefan-Boltzmann limit in the free theory has a known 
dependence on N t : e SB = 3p S B = (N* - 1)7T 2 /15. R(N t ) where R(N t ) = 1 + Sn 2 /21N? + ... is the 
known discretization error in the integral method [ |l4"[ HQ. When comparing data at different N t in 
Fig. ||, we normalize them by this known discretization error. 

As Fig. U indicates, there is a considerable cutoff effect in e — 3p in the temperature range 
< 1 . 1 T c for the N t = 6 lattice, for both SU(4) and SU(6). For T > 1 3T C for N, = 6, the cutoff effect 
is small, which is also expected since a = 1 /l.ST c and is near the scaling regime. For pressure, the 
cutoff effect is already seen to be small at T c for N t = 6, for both the theories. We also see that 
within the statistical uncertainties, volume dependence is small already for an aspect ratio £ = 3. 

In Fig. ^ we show the results for SU(A^ f ) thermodynamics for N c = 3, 4 and 6. For Af c = 3, we 
have taken the plaquette data from Ref. Jl3| ] and analyzed it similarly to N c = 4 and 6. To emphasize 
both the deviation from the Stefan-Boltzmann limit, and to look at evidence of corrections to the 
leading Af? behavior of the extensive quantities, we plot the energy density and pressure normalized 
by their Stefan-Boltzmann values, e/esb and p/psb- Two comments are in order here. First, 
The thermodynamic quantities differ from their Stefan-Boltzmann values considerably even for 
temperatures near 4 T c . Also, the approach to the Stefan-Boltzmann limit is slow. Second, the 
correction to the leading N% behavior is small, as there is no systematic change between N c = 3 to 
6. Note that a subleading correction is expected to be /N*), and therefore, four times larger in 
SU(3) than in SU(6). 
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Figure 4: Results for energy density and pressure for SU(3-6) gauge theories, normalized to the correspond- 



ing Stefan-Boltzmann values. For SU(3), we have used the plaquette data of Ref. [ 13 1 



4. Summary and Discussion 

In this report we have studied the deconfinement transition in S\J(N C ) gauge theories, with 
emphasis on the continuum results. We have studied in detail the transition point, using a finite 
volume analysis to establish the first order nature of the transition for N c = 4 and 6. We have also 
found that for > 8, the two-loop RGE works quite well at the transition point. We have used 
this to set the scale of the theory, by transmuting the coupling at the deconfinement transition, j8 c , 
to get T c /A„s. We find this quantity to be very weakly dependent on N c . 

Next we have studied the thermodynamics of the deconfined plasma. Our results for scaling 
imply that for calculating the thermodynamic quantities near T c , one needs lattices with A^ > 8. It 
was found that an aspect ratio of 3 was enough for the finite volume effects to stay under control. 

Our estimates for the continuum thermodynamic quantities are shown in Fig. ^. The figure 
denotes that the correction to the leading Af? contribution is rather small in both energy density and 
pressure, already for N c = 3. It is also found that the energy density and pressure differ significantly 
from their Stefan-Boltzmann values, even deep in the deconfined phase. 

The large deviation from the Stefan-Boltzmann value, and the fact that some strongly coupled 
conformal field theories show similar deviations from Stefan-Boltzmann limit [0], have sometimes 
been used in the literature to speculate about a strongly coupled, conformal regime in pure S\J(N C ) 
gauge theories. The thermodynamics results can be used to investigate the feasibility of such a 
phase. Following Ref. [p3[], in Fig. |we plot the energy density vs. pressure, normalized by the 
corresponding Stefan-Boltzmann values. By construction the point at (1,1) is the Stefan-Boltzmann 
limit, while the diagonal line denotes conformality. Also shown are the weak coupling lines for the 
theories with the different number of colors [|l6|]. We find that the weak coupling line is reached 
before the conformal line, indicating the absence of a strongly coupled, conformal phase in the 
S\J(N C ) gluon plasma. 

We thank Mikko Laine for providing us with the weak coupling results of Fig. |5[ The com- 
putations were carried out on the workstation farm of the department of theoretical physics, TIFR. 
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Figure 5: Approach to conformality in SU(N C ) gauge theories. Also shown are the weak coupling results 
from [|16[. 

We thank Ajay Salve for technical support. 
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